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■ Abstract. Time evolution of wave packets built from the eigenstates of the Dirac 

t— I | equation for a hydrogenic system is considered. We investigate the space and spin 

motion of wave packets which, in the non-relativistic limit, are stationary states with 
f~>) . a probability density distributed uniformly along the classical, elliptical orbit (elliptic 

^5 WP). We show that the precession of such a WP, due to relativistic corrections to the 

energy eigenvalues, is strongly correlated with the spin motion. We show also that the 
1 1 ^ , motion is universal for all hydrogenic systems with an arbitrary value of the atomic 

£h ' number Z. 
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PACS numbers: 03.65.P, 03.65, 03.65.S 
The detailed study of the time evolution of quantum wave packets (WPs) in simple 
atomic or molecular systems has been the object of growing attention for more than 
ten years [|]]. Most of the previous theoretical studies were done in non-relativistic 
framework. In the field of relativistic quantum mechanics most of the efforts have 
been focused on the problem of the interaction between the atoms and a mixture of 
static fields with, most of the time, intense laser fields @-ll|. Under these conditions 



the use of a relativistic theory is fully justified since the external field is then able to 
bring considerable energy to the WP. For isolated atoms, however, the use of relativistic 
dynamics is more questionable, if the WP is followed or observed only during a short 
period of time. In ref. Jl2| relativistic wave packets, corresponding to circular orbits, 



have been constructed for hydrogenic atoms with a large Z, and propagated over time 
according to the Dirac equation. Particular attention was paid to the spin collapse 
event, i.e. to the maximum entanglement, in the course of time, of the spin degree of 
freedom with the spatial ones. This phenomenon was first shown to take place for a 
WP in a harmonic oscillator with a spin-orbit force fI3|| , where it is periodic. For this 
reason it has been called the spin-orbit pendulum. In the Dirac equation with a Coulomb 
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potential, it is produced by the built-in spin-orbit force and it is not periodic. The time 
scale where this effect can manifest itself was discussed in ref. ]nj, as a function of 
the charge Z of the atom and the average principal quantum number N of the WP. 
This phenomenon is expected to take place on a longer time scale like the other time 



dependent quantum effects of spreading and of fractional revivals HT3fl . We intend below 
to complement this work by studying the possible relativistic precession of elliptic WPs 
and by comparing this precession to the spin motion. A preliminary version of the 
present work was already presented at the ECAMP VII conference [15[ . 



There are two possible ways to build up a WP 'on top of a classical elliptic orbit' 
in non-relativistic mechanics. One of them is by kicking properly a well designed WP, 
for example a Gaussian WP, that is then evolved in time by the free Hamiltonian until 
it spreads above an average classical ellipsis. This method is not very easy to apply and 
its main inconvenience is to produce an internal motion within the ellipsis that is able 
to blur the precession. Therefore, we have preferred a second method, much simpler 
and even more elegant, which consists of using the coherent WPs of ref. [|16|], which are 



stationary states of the non-relativistic Coulomb problem. The space probability density 
of these states was indeed shown to be highly concentrated around a classical Bohr- 
Sommerfeld ellipsis. If these states can be adapted to relativistic dynamics their time 
evolution will doubtlessly be due to relativistic effects, i.e. due to the fine structure that 
will act as a perturbation. Let us first show how to adapt these states to a relativistic 
theory. We want the spatial part of the large components of the new WP to tend (in 
the non-relativistic limit) toward the state \nj) of JTBJ which is defined as 

(Z + m)!(Z-m)!(2Z + l)l 1/2 



{n - I - \)\{n + l)\ 



I \ = y-/_ n a+m)/2 2 " 1 l( y n l ! 

' V U { ' [(l-m)/2]\[(l + m)/2}\ 

V 2/ V C ° S 2y \n,l,m) = \n,l,m) . (1) 

The probability density (717I717) is fairly localized onto a Bohr orbit with 
eccentricity e = sin 7 and the average angular momentum is 

Lv = ("-!) cos 7 , (2) 

where n is the principal quantum number of the orbitals \n, I, m) which are admixed 
in ([!]). The sum contains n(n + l)/2 values of m with I + m even. The contribution 
of states with m < I decreases very rapidly with m (more than one order of magnitude 
for each 2 units of m as shown in fig. 1). The dominant weights are those with 
m = I and their distribution is nearly Gaussian. The relative contribution of states 
with m < I increases, however, for larger values of the eccentricity parameter e. The 
larger admixture of these states causes much faster daphasing of the WP. Therefore for 
illustration of typical precession phenomena he have chosen the case e = 0.4. 

In order to study the entanglement of the spin degrees of freedom with the orbital 
ones, it is natural in a non-relativistic theory to start from a product state of an arbitrary 
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spinor y ^ j with the state \nj). 

I*«r> = |"7> (°)=£w£?M,™> ( j) ■ (3) 

It may be expanded in eigenstates of total angular momentum \n, I, j = I + s, rrij) with 
rrij = m + 1/2 or m — 1/2 and s = +1/2 or —1/2. 



|*nr> = W lm < " 



Im 



+ b 



>l + 1 + m i ix I — in , 

KZ,j>,m + 1/2)- W— — KZ,j<,m + 1/2) (4) 



2/ + 1 



21 + 1 



// + 1 — m l + m. , . , 

|w,/,J>,m- 1/2) + y — : - \n,l,j <1 m- 1/2) 



2/ + 1 



2Z + 1 



In a similar way as in ref. ||12|| , for a circular WP, the state (Q) is transformed into 
a four component relativistic state \l/ r by replacing in @ the non-relativistic states 
\n,l,j,mj) by the eigenstates of the Dirac equation with the same quantum numbers. 
In this manner the WP gets small components in the most natural way. The radial 
parts of the large and of the small components are taken obviously as different ones. 
The probability density of the relativistic WP built in this way is represented in fig. 2. 

The time evolution of the WP is produced by introducing the phase factors 
exp (—iEi~ t) and exp (—iEf t) as coefficients of states with j — 1 + 1/2 and j — I — 1/2 
with their corresponding eigenvalues. The four components Cj(t) i — 1, . . . , 4 of the WP 
at time t are given in the [Appendix A| with \l/ r defined as 

/ [ci(*)> \ 
|c 2 (t)> 

Mt)> 

V |c*(f)> J 

It should be stressed that \l/ r for t = is not any more a product state of the form of 
eq. (0) due to the built-in entanglement contained in the solutions of the Dirac equation. 
However, since the small components of \I/ r are indeed very small, this defect has no 
important effect on the magnitude of the initial spin. Since the spin-orbit coupling effect, 
i.e. spin-orbit pendulum ]TJ , manifests itself more efficiently if s and I are perpendicular 



\%.(t)) 



(5) 



to each other, we choose for our discussion: 



b = -P 



i.e. 



1 

2 



(6) 



Precession of quantum elliptical states in the Coulomb field has already been 



considered in ref. [17], starting also from the same coherent state as done here. However, 
the precession was studied only in non-relativistic quantum dynamics as a perturbation 
effect and no treatment of the spin was attempted. 
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Let us finally discuss the time units relevant to our problem. To the lowest order 
in terms of the fine structure constant, the energy of an eigenstate \n, l,j = I + s) of the 
Dirac equation can be written (in a.u.) for a hydrogenic atom of charge Z as 

f -e 1 z4q2 (7) 

nl3 ~ n 2n\l + s + l/2) [) 

with s = +1/2 or —1/2. The constant energy E n produces no effect on WP, since it 
depends only on n. Let us define an average time unit T p (p for precession): 

2tt 

T P ~ / dE nll X ( 8 ) 

= 21 ™ Tk . (10) 

Brackets () in ©-© denote average values, l av is given by (§). For n = 50, e = 0.4 
precession time T p ranges from 1.96 • 10 _11 s for Z = 92 to 1.4 • 10 _3 s for Z = 1. Tk in 
(0) denotes the Kepler period. 

It is necessary to compare T p to the radiative lifetime T™f of hydrogenic levels. We 



will use the estimation of T^f for a single n, / state given in fL8| . In atomic units it is 



r//' 1 (/ • -) . (ii; 



2a 3 Z 4 V 2, 

which was found to agree within 10% with experimental data. One obtains the universal 
ratio: 

^- = ^^0.061 (12) 

which guarantees that the precession of the wave packet takes place long enough 
time before even a single photon is emitted. The occurrence of a in this ratio is 
understandable, since T p is a classical-like quantity, while T™f is proportional to 1/h, 
because it can be expressed as a ratio of a typical energy of the emitted photon to the 
classical radiation rate corresponding to the orbital motion. 

When t = T p the linear terms in the expansion of E n ij in the power of / contribute 
on the average to 2tt in the exponential factors and the WP is expected to restore. 
However the terms of the higher order dephase differently the various partial waves, 
and this leads to a spreading of the WP near the initial shape fl4| . See the discussion 
on the magnitude of these terms in |Appendix B| . Expression (10) or (11) (with the 
Kepler period TV = 2nn 3 /Z 2 ) is also recognized as the classical precession time in the 
relativistic Coulomb problem fl9|| . 

Let us note that dE n ij/dl is also (to the first order) the energy difference between 
two spin-orbit partners. Therefore the precession time can also be interpreted as the 
recurrence time of the spin. Hence we should observe a strong correlation between the 
spatial motion of the density: the precession, and the spin motion. 
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Let us discuss now the dynamics dependence on the atomic number, Z. Formula (|9|) 
suggests that the relativistic effects under consideration in this article, depend crucially 
on Z. For sure the highest possible Z are indeed required to lower T p as much as 
possible. It is interesting to stress nevertheless, that within a very good approximation, 
a scaling of Z is possible which leads to the universal behaviour of the wave packet ([[]). 

First of all since E, approximated by eq. (|7]), scales as Z 4 , i.e. as T~ x , the variable 
Z disappears from Et if we use the reduced time t/T p . The autocorrelation function 
( |A.20| ) which is expressed only in terms of Et is the simplest quantity which has a 
universal behaviour, provided the same n and e are used for all values of Z. 

The other quantities, like the probability density and the the spin expectation values 
depend on the radial wave functions and radial integrals. In a non-relativistic theory 
scaling of the radial wave function is elementary, it is obtained by dividing the radial 
wave fuction by Z 3 ^ 2 and multiplying the radial variable by Z. This is an exact property. 



The ratio small/large components of the relativistic solutions scale as J (1 — E)/(l + E), 
i.e. roughly like Z 2 and the radial wave functions depend also in a more complicated 
manner on Z. Nevertheless, on the whole as seen in fig. 5 below, the small components 
contribute a very small part of the probability density even for Z = 92 (see also fig. 1 
of ref . Jl2 ) . The scaling of the probability density is therefore entirely governed by the 
large components i. e. by the non-relativistic theory. 

In a similar way we have checked that the radial integrals which contribute to the 
spin expectation values have the same properties: the integrals G + ,G^, and G + _ are 
equal to 1, within less than 10~ 5 and the other F+,F- . . . contribute in a very small 
manner, also of the order of 10~ 5 . 

Therefore, the universal behaviour of the wave packet is justified and except for 
fig. 2 no value of Z is given. For longer times, the energy factors omitted in eq. (|7]) 
which involve higher powers of Z, play a role and produce a genuine Z dependence. 
Those effects will not be discussed here. 

The probability density of the wave packet with n = 50, e = 0.4 and a = b = l/y/2, 
with Z = 92, is shown for a set of times up to t = T p in fig. 3 and fig. 4. The part of the 
density coming from the small components, shown in fig. 5, is also concentrated on top 
of an ellipsis but its magnitude is a thousands times smaller than the total density. For 
t < T p / 4 the density precesses as described classically with a small dispersion. However 
the spreading takes place very rapidly for larger t and is quite extended for t = T p . 

The precession of the ellipsis, the recurrence and spreading can be seen in a more 
condensed manner in the autocorrelation function represented in fig. 6 for three different 
spin preparations (spin up, spin down and a = b = 1/^/2). For small t/T p the WP 
becomes almost orthogonal to its initial parent, for t = T p a recurrence occurs but the 
overlap is only 0.7. Another peak occurs at t = 2T p but higher frequencies become 
important and spread the recurrence. For t > 3T P these higher frequencies play a 
dominant role. An example of an approximate revival for t = 22.6 T p is displayed in 
fig. 7. Fractional revivals [I4"| can also be seen to some extent. Two examples of 1/3 
and 1/2 revivals are presented in fig. 8. 
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The rough independence of the autocorrelation function on the spin preparation 
requires some explanation. We can approximate this function by 

(* r (0)|* r (t)> « a 2 ]>>2 exp (-iE+t) + ^E^exp {-iEft) , (13) 



i i 

,2 



where we have neglected the terms with very small weights wf m ^i. The contribution of 
the two sums above is almost the same because wu has a smooth variation with I on the 
one hand, and because of the exact equality E[~ = on the other hand. 

Finally the time evolution of the spin averages is presented in fig. 9. This figure 
exhibits what we can call the relativistic spin-orbit pendulum. Although the wave 
packet is not prepared initially in a pure state of spin, the impurity is very small (for 
a = b = l/v2 one has (a x ) ~ 0.9997). As time goes on the spin stays very nearly in 
the Oxy plane and rotates around Oz with period T p . Its magnitude slowly decreases, 
however, and for 5 < t/T p < 10 the spin is amost totally entangled with the orbital 
motion, since the average of its three projections are almost zero. During a period of 
time that last for about 5T p the angular momentum of the spin is transferred to the 
orbital motion and therefore the mean trajectory is not planar anymore (jOJ. Since the 
orbital angular momentum of 40h is much higher than h/2 this geometrical effect can 
hardly be represented graphically. From t ~ 15T P a revival of spin occurs. The spin 
rotates and increases its magnitude at the same time. This event stays even longer than 
initially. However the recurrence is only partial. 

In conclusion we can say, that for not too long time, the precession and spin motion 
of the WPs are fairly well described by the following approximation: non-relativistic 
wave functions of the form (^j and relativistic energy eigenvalues. From this point of 
view the full, computationally very demanding, relativistic approach is unnecessary. 
However, this conclusion may be formulated not a priori but only a posteriori. 

We would like to stress the richness of the dynamics just described. Indeed, in 
addition to the relativistic precession of the ellipsis we have obtained for longer times 
its fractional revivals. During the evolution the spin of the electron is entangled with its 



orbital motion to various degrees. All this agree completely with previous results [14 



and [13|. However, it was obtained here in a full relativistic calculation. Since we have 
been able to scale the atomic number Z we have given a universal behaviour to our WP. 
It is, however, clear that this scaling is destroyed in real atoms in a more realistic theory 
which would take into account quantum defects. Their inclusion would also distort 
the dynamics, for example it would change the precession time, in a way that is out 
of reach of our simple theory. As far as purely relativistic effects are concerned, like 
the importance of the small components or the zitterbewegung, we have found them 
negligible in the Coulomb problem, in contrast to the Dirac oscillator in which they 
play a major role. 
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Appendix A. Details of the calculations 

Replacing a non-relativistic wave function in (|4]) by the eigenstates of the Dirac equation, 
one obtains for ^ r (t): 



q (t) - V w {n) aJ l + 1 + m ( ^'^1,^+1/2 



tin 



+a \hr-^\ t o exp(-iE,t) (A.l) 



2/ + 1 \ — /n'_^i-l,i<,m+l/2 

l l + l- m ( ig n > + £li t j >>m -i/ 2 

21 + 1 \ —f n ' + &'l+l,j > ,m-l/2 
11 + 




V 2/ + 1 y —fn'_Ql-l,j < ,m-l/2 , 

We have used the notations of [|12|: g(r) and /(r) are the radial parts of the large and 
small components associated with the quantum numbers n' = n — (j + 1/2). These 
functions are multiplied by the spherical tensors Fli,j, m which are defined by eq. 4a and 



4b of ||12|| . The energy of the spin orbit partners of a given value of I is denoted by 
if j = I + 1/2 and Ef if j = I — 1/2, respectively. 

For numerical calculations it is convenient to rewrite components of (|A.1|) in the 



following form 

/ \ \ I / i-i 4- \ (n) I I + 1 + TYl . ^ 

ciOO) = 1 2^ yn> + exp {-lEft) ^ w K J - - - \l, m) 



J(l + l-m)(l + m) 
+b^ -|/,m-l> I (A.2) 



+ g n >_ exp {-iE l t) w\m ( a ir7T^ m ) 



m 



21 + 1 

J(l + l-m)(l + m) 
~ b ~ 2p~* \l,m-l) 



|ca(*)> = i J2 {#< exp (-iE+t) £ U— ti-^|/, m ) 



i/(Z + 1 + m)(/ - m) 
+a i^ ^ -|Z,m + l>| (A.3) 



21 
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+ q n ,_ exp {-iE l t) J2 Am (^T^J % m ) 

yj(l + 1 +m)(l - m) 



21 + 1 



\l,m+l) 



|c 3 (t)) = EU«v ex p(-^)E^t ) 



(l + l + m)(l + l-m) 
\J (2/ + l)(2/ + 3) 



\l + l,m) 



+b A 



(l + l-m)(l + 2-m) 
\ (2/ + l)(2/ + 3) 



|Z + l,m-l) (A.4) 



+ f n >_ exp (-iEi t) ]T 



(n) 



(/ + m)(Z — m) 
\ (2Z + 1)(2Z- 1) 



|Z - l,m) 



-6* 



(Z + m)(Z -1 + m) 
\\ (2/ + l)(2/-l) 



|Z- 1,771- 1) 



M0> = E ex p (-^*) E «ff ( 



(Z + 1 + ra)(Z + 1 - m) 
^ (2Z + l)(2/ + 3) 



Z + l,m) 



(Z+ 1 + m)(/ + 2 + m) 
\| (2Z + l)(2Z + 3) 



|Z + l,m + l) (A.5) 



+ /„/_ exp {-iE l t) E 



(n) 
Im 



(I + m)(Z — m) 
\ (2Z + 1)(2Z- 1) 



|Z - l,m) 



(Z — m)(l — 1 — m) 
\| (2Z+1)(2Z-1) 



Z - 1,771 + 1) 



Using the explicit form of the WP, one obtains for the average values of the spin 
operators: 



(cr x )t = 2a&E 



w 



l,m 



Im 



(Z + l) 2 -m 2 ri l 2 -m 2 r (Z + l) 2 -m 2 



-F. 



(2Z + 1) 2 
l 2 -m 2 



(21 + 1)(2Z- 1) 



2G 



(21+ I) 2 ^(21+ l)(2Z + 3) 
Z(Z + l)+m 2 



(2Z + 1) 2 



cos (uit) 



+Wl m Wl- 2 ,m-2 



(Z + m)(Z - 1 + m)(l - 2 + m)(l - 3 + m) 



(21- 1) 2 (2Z + l)(2Z + 3) 

y/(l + m)(l -1 + m) 
(2Z + 1) 



cos (ou'/t) 



(A.6) 
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(l + 2-m)(l-3-m) y/(l + 2-m)(l + l-m) 

-T — — ; — : ; 1 a. 



(21 - 1) 



(21 + 3) 



+Wi !m W l+2 ,m—2 



2F_ 



(/ 2 -m 2 )((/-l) 2 -m 2 ) 
\l (2Z-l) 2 (2Z + l)(2Z-3) 



cos (w'l't) 



F' 



1 



(Z + 1 - m)(l + 2 - m)(l + 3 - m)(l + 4 - m) 



(21+ l)(2Z + 3) 2 (2Z + 5) 



cos (u^'i) 



lm l 

-Wl,mWl-2,m-2 

+ W, m W[_2 m 2F_ 



2m 
2/ + 1 

F. 



G + _ sin (u^t) 



(A.7) 



A 



(I + m)(l - 1 + m)(l - 2 + m)(l - 3 + m) 



(21- 1) 2 (2Z + l)(2/ + 3) 



sin (uj'/t) 



(I 2 - m 2 ) ((I - If - m 2 ) 
\J (2Z- 1) 2 (2Z+ l)(2Z-3) 



sin (uj'(t) 



-Wl :m Wl +2 , m -2 

and 



F' 



A 



(Z + 1 + m)(Z + 2 - m)(Z + 3 - m)(l + 4 - m) 
(21- l)(2/ + 3) 2 (2/ + 5) 



sin (u^'i) 



2 2m + 1 / Z + 1 +m 
a — — — I G + — ^ — : G 



2Z + 1 



+b 2 ^-^ I G_ 



2/ + 1 
/ + 1 — m 



I — m I + 1 + m I - 

2l + l~ + 21 + 3 + ~2T 



2Z + 1 



2/ + 1 



G. 



I + m 
2Z + 1 



F + l+ J~ m +K l + 



2Z + 3 



2Z 



m N 



cos (uj(t) I a z 



(Z-m)(Z + l + m) 2 (/ + m)(Z + l - m)' 



(21 + If 



-b 



+4(a 2 - b 2 ) w hm w l+2 , m F'_. 



((I + l) 2 - m 2 )((l + 2) 2 - m 2 ) 



■\ (2Z + l)(2Z + 3) 2 (2Z + 5) 

In the above formulas, the following notations have been introduced: 

^ = (E+ - Ef) , 

uj[ = (Er +2 - E+) , 

J{ = (Ef - BU) 
Note that uj[ = uo'i +2 - Radial integrals are denoted as follows: 

poo . . 2 

G + = I (g+(r)) r 2 dr , 

poo , . 2 

G- = j o (gf(r)) r 2 dr , 
F + = jH(/, + (r))Vdr, 



(21 + If 
cos (oj[t) 



(A.8) 

(A.9) 
(A.10) 
(A.H) 

(A.12) 
(A.13) 
(A.14) 
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gii. r )gii. r ) r2 dr , 

i 

fit 2 {r)fr{r)r 2 dr , 

i 
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(A.16) 
(A.17) 
(A.18) 
(A.19) 



Apart from the case of G + , GL, F + , _F_ for I = n — 1, which are relatively easily obtained 
analytically, all other radial integrals have been calculated numerically (using quadruple 
precision). 

The autocorrelation function can be calculated from (|A.1|) in a straightforward way: 



(*V(0)|¥ r (t)) = £ expHE+t) 



+ exp {—iE l t) 



E< 



2 Z + 1 + m 2 / + 1 — m 



21 + 1 



2/ + 1 



E 



2 I 2 l — m A + m 
Wl ' m I 2TTT + 2/TT 



(A.20) 



Appendix B. 

Let us use the following notation: 
k = j + l/2 , S k 



E 



nlj 



m c z 



X 



(Za) 



The exact eigenergies £k (in units m c 2 ) are given by 



1 + 



X 



In 



-k + ^/k 2 ^ 



X 



2\2 



-1/2 



Expanding this expression in Taylor series with respect to x one obtains 



1 - 



X 

2rt 



X 



—(- 



3 

2n 



x 



+ 



4n 3 \2k 3 2nk 2 n 2 k An 3 



3 +A) + [0(x)] 1 



(B.l) 



(B-2) 



(B.3) 



In eq. |7| only the lowest term depending on k in x A i.e. 5^8^ = —fsj: has been 
included. The higher order term 8q8^ contributes very little, because the ratio 
5Q8k/5^8k = x 2 (-^2 + 4^ — 2^2 + jpg-) reaches the maximum value about 0.0005 for 
e = 0.4 and Z = 92 and stays much smaller for lower Z. Then the time evolution for 
not too long period is mainly determined by the lowest order contribution fl7|). 
The precession time is determined by the derivative 



d£ k 



x 



dk 



2n 3 k 2 



, 3x 2 ( 1 
1 + I - 



(B.4) 



2 \2k 2 nk 2n 2 J _ 

Again in eqs. only term of the order of x A has been used. The x 6 -order term 

contributes at most about 0.00022 of the x 4 -order term for e = 0.4 and Z = 92. 
Therefore we conclude that the x 6 -order term can be safely neglected in estimation 
of the precession time T p . 
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Figure 9. Time evolution of expectation values of spin components. 



